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Members of the TMEM16 (Anoctamin) family of membrane proteins have been shown to be essential constitu- 
ents of the Ca 2+ -activated Cl~ channel (CaCC) in many cell types. In this study, we have investigated the electro- 
physiological properties of mouse TMEM16F. Heterologous expression of TMEM16F in HEK293 cells resulted 
in plasma membrane localization and an outwardly rectifying Ici,ca that was activated with a delay of several min- 
utes. Furthermore, a significant Na + current was activated, and the two permeabilities were correlated according 
to Pn 3 = 0.3 Pq- The current showed an EC50 of 100 uM intracellular free Ca 2+ concentration and an Eisenman 
type 1 anion selectivity sequence of P S cn > Pi > Pbi > Pa > Pas P - The mTMEM16F-associated Ici,ca was abolished in 
one mutant of the putative pore region (R592E) but retained in two other mutants (K616E and R636E). The 
mutant K616E had a lower relative permeability to iodide, and the mutant R636E had an altered anion selectivity 
sequence (Pscn = Pi = Pbi = Pci > Pas P )- Our data provide evidence that TMEM16F constitutes a Ca 2+ -activated anion 
channel or a pore-forming subunit of an anion channel with properties distinct from TMEM16A. 



INTRODUCTION 

Anion channels are involved in several physiological pro- 
cesses, including housekeeping functions such as mainte- 
nance of cytoplasmic ion composition, pH regulation, 
and cell volume regulation (Nilius and Droogmans, 2003; 
Hoffmann et al., 2009; Verkman and Galietta, 2009) and 
more specialized physiological funcdons such as trans- 
epithelial electrolyte/fluid transport in absorbing and 
excreting epithelia (Larsen, 2011) and muscle contrac- 
tion and excitability (Jentsch etal., 2002; Hartzell etal., 
2005). Furthermore, intracellular CP channels are im- 
portant for maintaining electroneutrality when protons 
and Ca 2+ ions are transported into intracellular compart- 
ments (Jentsch etal., 2002). These functions are main- 
tained by various types of CI channels whose mechanisms 
of activation and molecular identities are often not known. 
Activation stimuli include ligand-gating, cAMP-dependent 
phosphorylation, voltage-gating, cell swelling, and in- 
creases in intracellular Ca 2+ (Jentsch etal., 2002; Nilius 
and Droogmans, 2003; Eggermont, 2004). Ca 2+ -activated 
CI channels (CaCCs) were first described in Xenopus 
laevis oocytes (Miledi etal., 1982; Barish, 1983) and Sal- 
amander rods (Bader et al., 1982) and have since been 
found in many cell types where they have important func- 
tions in, e.g., membrane excitability in cardiac muscle and 
neurons (Andre et al., 2003; Guo et al., 2008) , olfactory 
transduction (Matthews and Reisert, 2003), epithelial 
secretion (Kunzelmann etal., 2007; Rock etal., 2009), 
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regulation of vascular tone (Angermann etal., 2006), and 
photoreception (Lalonde etal., 2008). 

At low Ca 2+ concentrations, CaCCs are characterized by 
time-dependent outward rectification and activation in a 
voltage-dependent manner by raised intracellular con- 
centration of free Ca 2+ ( [Ca 2+ ];) . At subsaturating [Ca 2+ ]i, 
CaCCs are slowly activated and deactivated by depolariz- 
ing and hyperpolarizing plasma membrane potentials, 
respectively. A higher Ca 2+ sensitivity at depolarizing po- 
tentials contributes to the outward rectification (Nilius 
et al., 1997). At saturating [Ca 2+ ];, the channels are fully 
activated at all physiologically relevant membrane poten- 
tials, and rectification is decreased with the current-voltage 
relationship approaching linearity (Kuruma and Hartzell, 
2000; Nilius and Droogmans, 2003; Hartzell et al., 2005). 
The [Ca 2+ ]; required for half maximal activation was 
reported to be ^200 nM (dependent on voltage) with a 
Hill coefficient of Ca 2+ binding > 1 (Nilius et al., 1997; 
Pedersen et al., 1998b; Kuruma and Hartzell, 2000; 
Klausen et al., 2007). CaCCs show an Eisenman type 1 
permeation profile given by Pscn > Pi > Pbf > Pci > Pas P , 
a single channel conductance of 0.5-5 pS, and sensitiv- 
ity to nonselective CI channel blockers like DIDS, 
NPPB, NPA, and niflumic acid (Jentsch etal., 2002; 
Nilius and Droogmans, 2003). CaCCs in some cells are 
activated directly by Ca 2+ , whereas CaCCs in other cells 
are dependent on the activity of CAMKII (Hartzell 
etal., 2005). 
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Recently, TMEM16A, officially named anoctamin 1 
(ANOl), was shown, by three independent laboratories 
to be associated with a Ca 2+ -activated CP current (Ici,ca! 
Caputo etal., 2008; Schroeder etal., 2008; Yang etal., 
2008). TMEM16A expression produced robust I C i,c a with 
Ca 2+ sensitivity, current-voltage relationship, ion selec- 
tivity, and pharmacology closely resembling canonical 
CaCCs. No obvious Ca 2+ -binding sites have been found 
in the protein (Hartzell etal., 2009). However, Yu et al. 
(2012) found that two amino acids, E702 and E705, are 
importantfor Ca 2+ binding. In the model of TMEM16A 
structure, these amino acids are placed in the extracellu- 
lar loop, and Yu et al. (2012) have therefore proposed a 
new structural model in which there is no reentrant loop 
and E702 and E705 are thus intracellular. TMEM16A 
belongs to the TMEM16 family of transmembrane 
(TM) proteins consisting of 10 members in mammals, 
TMEM16A-K (ANO 1-10). Hydropathy analysis pre- 
dicts members of the TMEM16 family to have eight TM 
domains, with the N and C termini located in the cytosol 
(Tsutsumi et al., 2005; Yang et al., 2008; Bolduc et al., 
2010; Milenkovic etal., 2010) The hydrophobic region 
between TM5 and TM6 seems to form a reentrant loop, 
which constitutes a part of the channel pore (Das etal., 
2008) , and mutations of positively charged amino acids in 
this region result in a reduction in the anion/ cation per- 
meability ratio (Caputo etal., 2008; Yang etal., 2008). It 
has been suggested that the TMEM 1 6s work as multimeric 
complexes (Schreiber et al., 2010), and recent data have 
shown that TMEM16A can form dimers (Sheridan et al., 
2011; Fallah etal., 2011), a functional structure it may 
share with other members of the family. 

Some TMEM 16 family members show tissue-specific 
expression in mammals. TMEM16B, TMEM16C, and 
TMEM16D are primarily expressed in neuronal tissue, 
and TMEM16E shows strong expression in skeletal mus- 
cle and the thyroid gland, whereas TMEM1 6F, TMEM 1 6H, 
and TMEM 16K have been found expressed in almost all 
tissue examined (Rock and Harfe, 2008; Stohr et al., 2009; 
Gritli-Linde et al., 2009; Schreiber et al., 2010) . TMEM16A 
shows broad expression but is most prominendy expressed 
in electrolyte-transporting tissue, neurons, and inter- 
stitial cells (Yang et al., 2008; Ousingsawat et al., 2009; 
Hwang et al., 2009; Schreiber et al., 2010). TMEM16A 
and TMEM16B have previously been shown to produce 
Ici.ca's, however, TMEM16B with somewhat different char- 
acteristics when compared with those of TMEM16A 
(Pifferi et al., 2009). TMEM16B is in particular expressed 
in photoreceptors (Stohr et al., 2009) and in the cilium of 
olfactory cells (Rasche etal., 2010), where it appears to 
play a role in sensory transduction (Pifferi et al., 2012). 

It is somewhat uncertain whether the remaining mem- 
bers of the TMEM16 family are CaCCs. Duran etal. 
(2012) recently published that TMEM16F and several 
other TMEM 16s are completely intracellular. In con- 
trast, Schreiber etal. (2010) found that expression of 



TMEM16F and TMEM16G in FRT cells induced a small 
ATP-stimulated I influx. However, the rate of I uptake 
by TMEM 16A was ^18-fold higher than the uptake rates 
by TMEM16F and TMEM16G, and a short isoform of 
TMEM16G, with intracellular localization and no pre- 
dicted TM domains, also took up I at roughly the same 
rate as the long isoform of TMEM16G and TMEM16F 
(Schreiber etal., 2010; Duran and Hartzell, 2011). In 
whole-cell patch clamp experiments, the authors reported 
a slowly activating Ca 2+ -sensitive membrane conductance 
in TMEM16F- and TMEM16K-expressing cells that was 
similar to that of mock-transfected FRT cells, within an 
8-min observation period (Schreiber et al., 2010) . Martins 
et al. (201 1) found that TMEM16F was an essential com- 
ponent of the outwardly rectifying CP current that can 
be stimulated by cAMP in the presence of CFTR. Re- 
cently, Shimizu etal. (2013) showed that TMEM 16F was 
a Ca 2+ -activated anion selective channel. In contrast, add- 
ing a great deal of confusion, it was recently reported that 
TMEM16F is a Ca 2+ -activated nonselective cation channel 
(Yang etal., 2012) and that TMEM16F is essential for 
Ca 2+ -dependent phospholipid scramblase activity involved 
in phosphatidylserine exposure to the surface of blood 
platelets. A patient with Scott syndrome caused by a de- 
fect in phospholipid scrambling activity had a mutation 
in TMEM16F (Suzuki etal., 2010; Castoldi etal., 2011; 
Nurden and Nurden, 2011; Yang etal., 2012). Finally, a 
study by Almaca et al. (2009) has suggested TMEM 16s to 
be implicated in cell volume regulation. However, this ef- 
fect seems predominantiy to be secondary to ATP release 
and purinergic signaling, resulting in CaCC activation 
(Kunzelmann etal., 2012), and Shimizu etal. (2013) 
found that neither knockdown nor overexpression of 
TMEM16F affected the volume-sensitive CP current in 
HEK293T (human embryonic kidney 293 T) cells. 

The present investigation aims at the function of 
TMEM16F in mediating a TM CP current. We show that 
mouse TMEM16F localizes to the surface membrane 
and generates Ici.ca when expressed in HEK293 cells. 
The current is activated after a relatively long delay (min- 
utes) in response to high [Ca 2+ ];. We determine the Ca 2+ 
dependence, relative anion/cation selectivity, anion se- 
lectivity sequence, and the effect of mutations in the pu- 
tative pore domain. 



MATERIALS AND METHODS 

Cell culture and transfection 

Two different model cell lines were used: the mouse-derived Ehrlich 
ascites tumor cells (EATCs), and the human-derived HEK293 
cells. The cells were kept at 37°C with 5% C0 2 /95% air. HEK293 
cells were kept in DMEM and EATCs in RPMI 1640 . 10% fetal bo- 
vine serum and 1% penicillin/streptomycin were added to each 
media. HEK293 cells were transfected by Trans-IT-293 (Mirus 
Bio LLC) and used for electrophysiological recordings 24-48 h 
after transfection. 
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Isolation of RNA, reverse transcription (RT), 
PCR, and cDNA cloning 

Total RNA was isolated from EATCs using the RNeasy Mini kit 
(QIAGEN) according to the manufacturer's protocol. cDNA was 
prepared in a total volume of 40 pi by hybridization of 500 ng oligo 
(dT) primers to 4 pg RNA at 65°C for 5 min, followed by extension 
at 42°C for 50 min in the presence of 200 U Superscript II RT (In- 
vitrogen), 500 pM dNTP, 10 pM DTT, 50 mM Tris-HCl, 75 mM KC1, 
and 3 mM MgCl 2 , pH 8.3. Finally, the RT was inactivated at 70°C for 
15 min. PCR was performed in a total volume of 20 pi containing 
1 pi of the RT reaction, 0.5 mM dNTPs, 0.5 pM of each primer, 2 mM 
MgCl 2 , and 2 U Taq polymerase in PCR buffer (conventional PCR) 
or 2 U Pfu Ultra II fusion HS polymerase (Agilent Technologies) 
in Pfu PCR buffer (for cloning). PCR consisted of a denaturizing 
step of 2 min at 95°C, followed by 35 cycles of 95°C for 30 s, 58°C 
for 30 s, and 72°C for 2 min. Primers used for PCR were 16A_for, 
5 '-GACCTGGGCTATGAGGTTC A-3 ' ; 16A_rev, 5 '-GGCTGATGTC- 
TTTGGGGATA-3 ' ; 16B_for, 5 '-CCAGGGAAGCAGAGTTCTTG-3 ' ; 
16B_rev, 5 '-TGTTGTTGGCTCGAGAACAG-3 ' ; 16C_for, 5'-TGA- 
CAAGTC AGCGTTTCCAG-3 ' ; 16C_rev, 5'-CAGCCAGGCGAAGTA- 
TAAGG3'; 16D_for, 5 '-TTTGAGCACCTCGTGTTTTG-S ' ; 16D_rev, 
5 '-TCATGGCCACTCATTGTGAT-3 ' ; 16E_for, 5'-TCTGGATCATG- 
CCTGAACTG-3' ; 16E_rev, 5'-ATGGTTAGCTGGGTGGTCAG-3'; 
16F_for, 5 '-GCAGCCCTTGGATCTTATCA-3 ' ; 16F_rev, 5'-TGCTG- 
TAGCTCAACGGTGTC-3 ' ; 1 6G_for, 5 '-TCTTCCTGGTGGGATGT- 
TTC-3'; 16G_rev, 5 '-CCACAGTGCCATGAACAAAG3 ' ; 16H_for, 
5 '-ACGC AGGCAGAGCTAGAGAG-3 ' ; 16H_rev, 5' ACACTGGCC- 
AATGAGGTAGG3 ' ; 16J_for, 5 ' -TGAAGGATGGGGTTTTTGAG-3 ' ; 
16J_rev, 5 '-TGCC ATGAAGATGGCAAATA-3 ' ; 16K_for, 5'-TCT- 
GAGTGGACCAGCCTTCT-3 ' ; 16Krev, 5 '-AGAAGAGTGAGGCGA- 
AGCAA-3'; ARP_for, 5 '-CGACCTGGAAGTCCAACTAC-3 ' ; and 
ARP-rev, 5'-ATCTGCATCTGCTTG-3' . Real-time quantitative PCR 
was performed in triplicates using an MX4000 Real-Time PCR 
system (Agilent Technologies) and SYBR green PCR master mix 
(Applied Biosystems) in a total volume of 20 pi containing 1 pi of 
the RT reaction, 200 nM of primers, and 10 pi of 2x master mix. 
Primers used are listed above ( 1 6F and ARP) . Target expression level 
was normalized to the reference gene level (rl8S), and the relative 
expression ratio was calculated using the equation 

j-, AC(t) tar g el (control treated) 

„ _ ^target 
ratio ~ ~ AC(t) ref (control treated) ' 
^ref 

where E target and E re r are the PCR amplification efficiencies for the 
target and reference gene (ARP), respectively, and AC(t) target and 
AC(t) re f are the change in C(t) values for target and reference 
genes. Primers used to clone and generate attB-flanked gateway- 
compatible PCR products were 16F_for, 5'-GGGGACAAGTTTGTA- 
CAAAA^AGO\^GGCITCAC(^TGCAGATGATGACTAGGAAGGTG3'; 
1 6F_rev, 5 ' -GGGGACC ACTTTGTACAAGAAAGCTGGGTCTCAT- 
TCGAGTTTTGGCCGCACG-3 ' ; 1 6F_revfusion, 5 '-GGGGACCAG 
TTTGTACAAGAA^GCTG<XiTCTTCGAGTT^ ' ; 

1 6A_for, 5 '-GGGGACAAGTTTGTACAAAAAAGCAGGCTTCAC- 
CATGAGGGTCCCCGAGAAGTAC-3 ' ; 16A_rev, 5 ' -GGGGACC AC- 
TTTGTA(3AAGAAAGCTGGGTCCTA(3AGCGCGTCCCQVTGGGA-3'; 
1 6A_rev_fusion, 5'-GGGGACCACTTTGTACAAGAAAGCTGGGT- 
CCAGCGCGTCCCCATGGTACT-3' ; 16K_for, 5 '-GGGGACAAGTT- 
TGTACAAAA\AGCAGGCTirACrATGAGAGT^ 
and 16K_rev, 5 ' -GGGGACCACTTTGTACAAGAAAGCTGGGTCT- 
CAGGTAGCTTCCTTCCCATCT-3 ' . PCR products were excised 
from 0.8% agarose gels and DNA extracted using an E.Z.N.A gel 
extraction kit (Omega Bio-tek). PCR products of mTMEM16F 
were then recombined into pDONR 221 using BP Clonase accord- 
ing to the manufacturer's instructions (Invitrogen) to generate 
mTMEM16F entry clones. After transformation into Escherichia 
coli and selection using 50 pg/ml kanamycin, DNA was purified 



and inserts in entry clones confirmed by full-length DNA se- 
quencing. Finally, using LR-clonase mix, the inserts from entry 
clones were recombined into pcDNA6.2/EmGFP-Bsd/V5-DEST 
and pcDNA.3.1-DEST47 to generate vectors either coexpressing 
mTMEM16F and EmGFP or an mTMEM 1 6F-GFP C-terminal fusion. 
Expression plasmids were produced in E. coli using 100 pg/ ml am- 
picillin as selection agent. 

Site-directed mutagenesis 

The mTMEM 16F R592E, K616E, R636E, and I342A mutations 
were constructed using a multisite-directed mutagenesis kit (Agi- 
lent Technologies). Primers were R592E, 5 '-TTCCTGGGCAAAT- 
ACGAGAGCGAAGAGTGTG-3 ' ; K616E, 5 '-ATCATGGGGGGA- 
GAGGC AATCTGGAAC AAC-3 ' ; R636E, 5' CCCATGGGTTATGA- 
ATCTAATTGGAGAGTATAAAAGAGTCTCGGG-3 ' ; and I342A, 
5'-GGTCTGTGATCCTGACGCCGGTGGCCAGATCCTG-3'. PCR was 
performed in a total volume of 25 pi containing 100 ng pcDNA6.2/ 
EmGFP-Bsd/V5-mTMEM16F-DEST or pcDNA-mTMEM16F- 
DEST47(GFP-fusions), 0.5 mM dNTPs, 0.5 pM of primer, 2.5 pi of 
lOx reaction buffer, 0.75 pi of quick solution, and 1 pi of multien- 
zyme blend (Agilent Technologies) . PCR consisted of a denaturiz- 
ing step of 1 min at 95°C, followed by 30 cycles of 95°C for 1 min, 
55°C for 1 min, and 65°C for 16 min. The PCR reaction was treated 
with Dpnl for 10 min at 37°C before being transformed into E. coli 
and selected for ampicillin resistance (100 pg/ml). Base pair sub- 
stitutions were confirmed by DNA sequencing. 

Electrophysiological recordings 

Cells were plated on poly-L-lysine-coated coverslips. Whole-cell 
voltage-clamp recordings were performed with the Axopatch 200B 
amplifier interfaced to a Digidata 1440A using pClamplO for re- 
cording and analysis (Molecular Devices). Analogue signals were 
acquired at 2.5 kHz and filtered at 1 kHz. All recordings were per- 
formed at room temperature (20°C). Patch pipettes were fabri- 
cated from borosilicate glass capillaries using a DMZ-Universal 
Puller (Zeitz Instruments) with a resistance of 2-3 Mfl when filled 
with the internal solution. For activation of Ici,c a > the pipette solu- 
tion contained (in mM): 100 Cs-Aspartate, 40 CsCl, 1 MgCl 2 , 1.5 
Na r ATP, 5 EGTA, and 10 HEPES, pH 7.2 with CsOH, osmolarity 
adjusted to 290 mosmol/kg. EGTA stock solutions were prepared 
by NaOH titration giving a final Na + concentration in the pipette 
solution of 15.3 ± 0.6 mM as measured by flame spectrophotometry 
using a FLM3 flame photometer (Radiometer). CaCl 2 was varied 
between experiments to investigate the role of 0, 0.19, 4.86, 4.98, 
5.11, 5.23, 5.41, and 5.70 mM Ca 2+ . The free [Ca 2+ ]i was measured 
using a Ca 2+ -selective electrode (Perfection; Mettler Toledo) in 
combination with a DrDAQ data logger and Pico Logger software 
(Pico Technology). Calibration curve was performed by varying 
[CaCl 2 ] in a calibration solution of (in mM) 100 Cs-Aspartate, 
40 CsCl 2 , and 1 0 HEPES. Because of the strong buffer capacity of 
the pipette solutions, we observed super-Nernstian potentials at 
concentrations <1 pM as the result of Ca 2+ depletion of the elec- 
trode at standard conditions (Sokalski etal., 1999). Thus, the free 
[Ca 2+ ] for low-level CaCl 2 was measured using the Ca 2+ -sensitive 
fluorescent probe, Fura-2P (Molecular Probes), as previously de- 
scribed (Pedersen et al., 1998a). The superfusate was a modified 
Ringer's solution and contained (in mM) 140 NMDG-C1, 10 NaCl, 
1.5 CaCl 2 , 1 MgCl 2 , 10 HEPES, and 10 glucose, pH 7.4 with NaOH, 
osmolarity adjusted to 315 mosmole/kg. In the anion selectivity ex- 
periment, NMDG-C1 was substituted by equimolar concentrations 
of NaCl, Nal, NaBr, NaSCN, or Na-aspartate (applied in a random- 
ized order) . The reference electrode was an agar bridge composed 
of 3% agar and 97% of the above mentioned superfusate solution 
containing NMDG-C1. Liquid junction potentials (Vy) were calcu- 
lated using the Clampex software (Molecular Devices) , and V test was 
corrected (V m = V test - V,j) in all the figures by - 19.3 mV, - 14.3 mV, 
-11.9mV, -14.4mV, -14.7mV,and 1.9mVfor NMDG-C1, NaCl, 
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NaSCN, Nal, NaBr, and NaAsp, respectively. Series resistance was 
compensated by 60-70%. A ramp protocol ranging from V te st = 
— 100 to 100 mV over 1 s with a holding potential of V test = 0 mV 
was run continuously with 15-s intervals to follow the current ac- 
tivation. A step protocol with steps ranging from V test = — 100 to 
100 mV with 20-mV increments and a step length of 1 or 5 s was 
recorded at peak I C i,ca- 

Immunofluorescence 

For images in Fig. 2, HEK293 cells were seeded on poly-L-lysine- 
coated coverslips in 6-well plates and grown until 50-60% conflu- 
ence. Cells in each well were transfected with mTMEM16F-GFP 
fusions using 5 pi PEI transfection reagent and 1 pg DNA. After 
24-h incubation, the transfection medium was replaced with fresh 
medium and the cells were incubated overnight. The cells were 
then washed in modified Krebs solution and incubated for 20 min 
at 37°C with modified Krebs solution containing 500 pM free Ca 2+ 
and wheat germ agglutinin (WGA; 1:200). The cells were washed 
in PBS, fixed with 2% paraformaldehyde in PBS for 15 min, washed 
twice in PBS, and permeabilized for 10 min with 0.1% Triton X-100 
and 1% BSA in PBS. The permeabilized cells were washed and 
then blocked in 2% BSA in PBS for 40 min and then labeled with 
primary antibody (1:300 rabbit anti-GFP in PBS with 2% BSA) for 
45 min. The cells were then washed three times in PBS with 2% 
BSA and incubated with Alexa Fluor 488 (1:300 in PBS with 2% 
BSA) for 45 min. After two washes in PBS with 2% BSA, the cells 
were incubated with DAPI (1:600 in PBS with 2% BSA) for 5 min 
and then washed three times in PBS. Finally, the slides were 
mounted (90% glycerol, 10% of lOx PBS, and 2% wt/vol JV-propyl 
gallate). Images were recorded using a BX-63 epifluorescence mi- 
croscope (Olympus). Primary antibody was GFP (A6455; 1:500 
Life Technologies). Secondary antibody was Alexa Fluor 488 
(A11055; 1:1,000; Life Technologies). For images in Fig. 4 E, 
mTMEM16F-GFP was expressed in HEK293 cells. Immediately be- 
fore fixation, the cells were exposed to ionomycin for 20 min in a 
solution containing 500 pM free Ca 2+ . The images were taken using 
an epifluorescence microscope. mTMEM16F-GFP was visualized 
with anti-GFP and Alexa Fluor 488 (green), whereas the plasma 
membrane was stained using rhodamine-conjugated WGA (red). 
Nucleus was stained with DAPI (blue) . Images in Fig. 6 F were re- 
corded in the EasyRatioPro software (PTI) using an Eclipse TI in- 
verted microscope (Nikon) equipped with CoolSNAP HQ2 digital 
camera (Photometries) and Deltaram X monochromator (PTI). 
HEK293 cells were seeded in 6-well plates and grown until 50-60% 
confluence. Cells in each well were transfected with mTMEM16F- 
GFP fusions using 6 pi Mirus293 transfection reagent and 2 pg DNA. 
After 24-h incubation, cells were trypsinized and 200 pi of cell sus- 
pensions transferred to new wells containing poly-L-lysine-coated 
coverslips. After 1 h, the settled cells were fixed in a final concen- 
tration of 2% paraformaldehyde, washed in PBS, and permeabi- 
lized for 15 min in 0.1% Triton X-100 in PBS. The permeabilized 
cells were washed, blocked in 2.5% BSA for 30 min, and labeled 
with primary antibody (1:300 rabbit anti-GFP in PBS with 1% 
BSA) overnight. The cells were then washed three times in PBS 
with 1% BSA and incubated with FITC-conjugated anti-rabbit sec- 
ondary antibody (1:400) for 2 h. After three washes in PBS, slides 
were mounted (90% glycerol, 10% of lOx PBS, and 2% wt/vol 
jV-propyl gallate) . 

Surface protein biotinylation 

Isolation of plasma membrane proteins used for the measure- 
ments in Fig. 6 G was performed using a Pierce Cell Surface Pro- 
tein Isolation kit (Thermo Fisher Scientific). In brief, cells where 
washed in ice-cold PBS and incubated for 30 min in PBS with 
250 pg/ml Sulpho-NHS-SS-Biotin. The reaction was quenched 
and cells were collected. After a single wash, cells were dissolved 
in lysis buffer and disrupted by sonication. 50 pi of the clarified 



lysate (10,000 g, 2 min, 4°C) was collected for protein concentra- 
tion measurement (DC protein assay [Bio-Rad Laboratories] using 
BSA as a standard) and total protein analysis. As protein concen- 
tration did not differentiate between samples, 600 pi of clarified ly- 
sate was incubated for 1 h with NeutrAvidin agarose beads in a 
centrifugation column on an end-over-end mixer for affinity puri- 
fication of biotinylated proteins. Beads were washed once, and 
protein elution was achieved using NuPAGE LDS sample buffer 
(Life Technologies) with 50 mM DTT Expression of specific pro- 
teins was analyzed using SDS-PAGE and immunostaining. 

SDS-PAGE and immunostaining 

SDS-PAGE was performed using the NuPAGE system (Life Tech- 
nologies) . Samples were diluted in LDS Sample Buffer (Life Tech- 
nologies) with DTT (total lysate = 42 mM; biotinylated membrane 
protein = 50 mM) . Proteins where separated by gel electrophoresis 
on a precast 10% Bis-Tris polyacrylamide gel (Life Technologies) 
using MOPS-SDS running buffer. Proteins were electrotransferred 
to Protran nitrocellulose membranes (Whatman), stained with Pon- 
ceau S red (Sigma-Aldrich). 



RESULTS 

Cloning of TMEM16F and localization experiments 
in HEK293 

TMEM16F was cloned from EATCs. As seen in Fig. 1 A, 
EATCs express mRNA for TMEM16A, TMEM16B, 
TMEM16F, TMEM16H, and TMEM16K In Fig. 1 B, this is 
also shown on the protein level for TMEM16A, TMEM16F, 
and TMEM16K by Western blotting with antibodies against 
A, F, and K. Full-length cDNA of the coding region of 
mouse TMEM16F was PCR amplified from EATCs. We 
also cloned the coding regions of mTMEM16A and 
mTMEM16K. As described in Materials and methods, 
cDNA was cloned using Invitrogen gateway technology 
and recombined into pcDNA6.2/EmGFP-Bsd/V5-DEST, 
which coexpresses EmGFP for easy detection of trans- 
fected cells, and into pcDNA3.1-DEST47 to generate a 
C-terminal GFP fusion for subcellular localization studies. 

Previous localization studies of TMEM16 proteins have 
reported plasma membrane expression of TMEM16A, 
TMEM16B, TMEM16E, TMEM16F, and TMEM16G (Bera 
et al., 2004: West et al., 2004; Rock et al., 2008; Schroeder 
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Figure 1 . TMEM16 expression in EATCs. (A) PCR for all known 
mouse TMEM16 genes in EATCs. (B) Western blots showing ex- 
pression of TMEM16F, TMEM16FA, and TMEM16FK protein in 
EATCs. Presented blots are representative of three blots made. 
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etal., 2008; Yang et al., 2008). In agreement with this, 
expression of GFP-tagged mTMEM16F in HEK293 cells 
showed significant staining at the plasma membrane 
(Fig. 2) . This is seen from an almost complete overlap be- 
tween the GFP-tagged TMEM16F and the plasma mem- 
brane marker WGA and no overlap with the ER marker 
Hsp47 (Fig. 2, right, merged images) . 

Whole-cell experiments of the l C i,c a 

in mTMEM16F-expressing HEK293 cells 

To study the putative Ici.ca's of mTMEM16F, we heterolo- 
gously expressed mTMEM16F in HEK293 cells and mea- 
sured whole-cell currents of cells showing expression of 
GFP. To suppress endogenous K + currents, the whole-cell 
voltage-clamp experiments were performed under con- 
ditions in which K + was replaced by Cs + in the pipette 
solutions. We did not observe activation of Ici.ca at 
submicromolar free [Ca 2+ ]i, butwith [Ca 2+ ]; in the range 
of 55-250 pM, a delayed and transiently activated cur- 
rent was identified. In the first set of experiments, 
mTMEM16F-transfected HEK293 cells were dialyzed 
with intracellular solutions containing a high free [Ca 2+ ]; 
(250 pM). Fig. 3 A (left) shows a representative current 
trace of the time course of the Ca 2+ -activated current de- 
termined by measuring the current at the beginning and 
at the end of the Ts voltage ramps from — 1 10 to 70 mV 
elicited at 15-s intervals (Fig. 3 A, left, inset). The outward 
current increased from 8.9 ± 1.2 to 86.1 ± 19.8 pA/pF 
(V m = 70 mV, n = 6, P < 0.05) within 7.9 ± 1.4 min. The 
voltage- and time-dependent current activated by 250 pM 
free [Ca 2+ ]; was measured at the peak of the current 
using 5-s voltage steps at potentials ranging from —120 
to 80 mV as shown in the right panel of Fig. 3 A Depolar- 
izing voltage steps (from 40 to 80 mV) resulted in out- 
ward currents with two components: an initial ohmic 
current followed by a delayed activation. Hyperpolariz- 
ing voltage steps (from —40 to —120 mV) also resulted 
in currents of two components: a small instant inward 
current followed by inactivation (not visible in Fig. 3 A 
because of scaling). Because of these time-dependent 



current changes, the steady-state currents become strongly 
outwardly rectifying with relatively small components 
athyperpolarizing potentials. Mock-transfected HEK293 
cells showed a similar but significantly smaller current re- 
sponse to whole-cell patch clamp with high [Ca 2+ ] but re- 
tained the same time-dependent activation/ deactivation 
at the peak of the response resulting in strong outward 
current rectification (Fig. 3 B). These currents are sup- 
posedly carried by the native human TMEM16F, which 
is expressed in HEK293 cells (Kunzelmann et al., 2009; 
Shimizu et al., 2013) . In our study, the variation in the re- 
sponse among the mock-transfected cells was fairly large, 
so the mean Ici.ca of a sample of five cells turned out to 
be nonsignificant with a change from 14.7 ± 5.9 to 26.6 ± 
5.0 pA/pF at 70 mV (V m = 70 mV, n = 5, P < 0.17) with the 
peak current at 6.3 ± 2.1 min. 

Importandy, we believe that the Ici.ca from mock-trans- 
fected cells shown here is not carried by an endogenously 
expressed TMEM16A, as expression of mTMEM16A gen- 
erated instandy activated, submicromolar Ca 2+ -sensitive 
currents similar to those reported in the original studies 
(Caputo etal., 2008; Schroeder etal., 2008; Yang etal., 
2008) . In fact, at an output gain of 1 mV/ pA and a [Ca 2+ ] ; 
of 1 pM, we lost control of TMEM16A-generated currents 
at potentials positive to 40 mV. The currents shown in 
Fig. 3 C were recorded at 70 nM [Ca 2+ ]i. 

Ca 2+ activation kinetics of mTMEM16F 

The Ca 2+ dependency of the mTMEM16F-associated Ici.ca 
was determined by measuring the above mentioned peak 
current as a function of [Ca 2+ ]iin the range of 55-415 pM. 
As shown in Fig. 4 A, with the typical time course of the 
current obtained by application of the step protocol, we 
observed a larger instantaneous current, a subsequent 
faster activation, and an increase in outward rectification 
of the CI current with increasing [Ca 2+ ]i- The resulting 
current-voltage relationships (Fig. 4 B) were obtained as 
described in Fig. 3 A. The mean Ici,ca densities at 70 mV 
were depicted as a function of [Ca 2+ ]i and fitted by the 
Hill equation giving an EC50 of 1 05.6 ±1.2 pM [Ca 2+ ] ; and 
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Figure 2. Subcellular localization ex- 
periments. mTMEM16F-GFP expressed 
in HEK293 cells. The images were taken 
using an epifluorescence microscope. 
mTMEM 1 6F-GFP was visualized with 
anti-GFP and Alexa Fluor 488. Anti- 
Hsp47 and Alexa Fluor 568 were used to 
label the ER, whereas the plasma mem- 
brane was stained using rhodamine-con- 
jugated WGA. The nucleus was stained 
with DAPI. 
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a Hill coefficient of 2.00 + 0.56 (Fig. 4 C). The CI" cur- 
rent was strongly outwardly rectifying both at low (90 uM) 
and high (190-250 uM) [Ca 2+ ]i and did not approach 
linearity at increasing [Ca 2+ ]; as seen for endogenous 
CaCC (Kuruma and Hartzell, 2000), TMEM16A, and 
TMEM16B currents (Yang et al., 2008; Pifferi et al., 2009; 
Scudieri et al., 2012). Rather, there was a tendency for the 
outward rectification to be more pronounced at increas- 
ing free Ca 2+ concentration (Fig. 4 D) , which was also ob- 
served for the endogenous Ici.ca- The currents showed 
strong time- and voltage-dependent activation at both 



low and high Ca concentrations (Fig. 4 A) . In spite of a 
fairly large random variation among the cells, it was indi- 
cated that the time from whole-cell breakthrough to acti- 
vation of the Ici.ca (compare with Fig. 3 A) depends on 
[Ca 2+ ] i, i.e. ; short delays with little variation were observed 
at the highest [Ca 2+ ]i (Fig. 4 D). To investigate the mech- 
anism of activation, we treated HEK293 cells expressing 
mTMEM16F-GFP with 10 pM ionomycin in an extra- 
cellular solution containing 500 pM Ca 2+ and made im- 
munofluorescence preparations after 20-min treatment 
(Fig. 4 E) . We did not observe any change in localization 
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Figure 3. Ia.ca densities in mock-, 
mTMEM16F-, and mTMEM16A-transfected 
HEK293 cells. (A) Typical current devel- 
opment recorded by a ramp protocol run 
every 15 s (left; see inset) and representa- 
tive current-voltage relationships (right) of 
mTMEM16F (250 uMfree [Ca 2+ ]i). Inward 
currents were initially inactivated similar to 
those in C; however, this inactivation is not 
visible at the given scale. (B and C) Mock 
(250 pM free [Ca^L; B) and mTMEM16A 
(0.07 pM free [Ca 2 *];; C). "a" indicates 
where the step protocols were recorded. 
(D) Mean peak I Ca ,ci current density when 
activated by 250 pM free [Ca 2+ ] ; of mock 
and mTMEM16F, measured at 70 mV and 
— 110 mV, respectively. Data are pre- 
sented as mean ± SE of 14 mock- and 8 
mTMEM16F-transfected cells and tested 
by one-way AN OVA: *, P < 0.05. 
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of mTMEM16F-GFP between nontreated and ionomy- 
cin-treated cells (compare with Fig. 2). 

Anion selectivity of mTMEM16F-induced currents 

The anion selectivity of mTMEM16F-induced Ici.ca was 
examined by measuring the relative permeability to vari- 
ous anions by replacing CI with equimolar concentra- 
tions of SCN~, Br , I - , or Asp - in the bath solution. The 



above bath solutions were prepared from Na salts of the 
respective anions. In a recent preliminary study, it was 
suggested that TMEM16F displays a cation conductance 
(Yang et al., 2011). In our study, the reversal potential of 
[Ca 2+ ];-activated membrane currents was significantly dis- 
placed toward positive values when NMDG-C1 of the bath 
solution was replaced by NaCl, which indicates a signifi- 
cant conductive Na + permeability (see below) . In a more 
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Figure 4. Ca 2+ activation kinetics of mTMEM16F. (A) Typical current-voltage relationships of mTMEM16F at max current with 90, 190, 
or 250 uM free [Ca 2+ ]j. (B) The mean mTMEM16F I-V relation at maximal activation with the indicated free [Ca 2+ ] i; generated from 
ramp protocol sweeps (from —120 to 80 mV). Data are presented as mean ± SE of three to nine replicates. (C) Hill plot of normalized 
mTMEM16F Ici.ca densities at 70 mV. Data are presented as mean ± SE of 3-12 replicates. (D) Ca 2+ dependence of mTMEM16F versus 
time from whole start of whole-cell recording to peak Ic a ,ci activation. Linear regression is included to illustrate the decline in the time 
to peak activation with increasing free [Ca 2+ ]i and is given by f(x) = — 0.04x + 14.11. (E) mTMEMlGF-GFP expressed in HEK293 cells. 
Immediately before fixation, the cells were exposed to ionomycin for 20 min in a solution containing 500 uM free Ca 2+ . The images were 
taken using an epifluorescence microscope. mTMEM16F-GFP was visualized with anti-GFP and Alexa Fluor 488, whereas the plasma 
membrane was stained using rhodamine-conjugated WGA. The nucleus was stained with DAPI. 
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detailed quantitative analysis, presented below, the ratio 
of the plasma membrane's Na + and CI permeability will 
be shown, P^/Pci ~ 0.3. Applying the Goldman-Hodgkin- 
Katz voltage equation for Pas P /Pci < 1> see below, and a 
PNa/Pci ratio different from 0, the expressions for the re- 
versal potentials before and after partial substitution of 
CI in bath with a monovalent anion, X~, are: 

V W - RT ln P cACn 0 +P Na [Na\ 
m -F n p cl [Cr] p +P Na [Na + l 

v m = RT Xn P c\Cr\+P x [X\+P Na [Na + \ 



PrACn t> +P N JNa + l 



Ris the universal gas constant (8.31 J/mol/K), Tis the 
absolute temperature (293 K), Fis the Faraday (96,485 
C/mol), P denotes permeability, subscript p refers to 
pipette solution, and o and s refer to bath solution be- 
fore and after substitution, respectively. With the mea- 
sured change in reversal potential given the symbol 



^■Vtw ~ ~ V**/' one proceeds as follows: 
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- + ln S- 



L [Na] 
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P ' P 



RT 



Exp I 



[CI 



>-[Na*] 




-FAV ra 
RT 
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'-[Na] 



[CI 



'-[Na] 



As [Na + ] p = 15.3 mM and P Na /Pci ~ 0.34 (see below), in 
Eq. 1, (PNa/Pa) x [Na + ] p ~ 5.2 mM, which amounts to 
<5% of [Cl~] 0 . Thus, as a first approximation, we ne- 
glect the term and obtain 



[cr\xE X p\^^\-[cr\ 



Fig. 5 A shows representative current-voltage relation- 
ships constructed from the linear voltage ramps of 
mTMEM16F-transfected HEK293 cells. Using Eq. 2, the 
shifts in reversal potential resulting from the change of 
anions in the bath solution were used to calculate the rel- 
ative permeability (P x /Pci; Table 1). The anion permea- 
bility sequence of the native Ici.ca in mock-expressing 
cells (putatively because of native hTMEM16F) was Pscn > 
Pi > PBr = Pci > Pas P (Fig. 5, A and C) , closely resembling 
an Eisenman type 1 selectivity sequence. The anion per- 
meability sequence of the Ici.ca in mTMEM16F-expressing 
cells was P SC n > Pi > Pbi > Pa > Pas P (Fig. 5, B and D) , an 



Eisenman type 1 selectivity sequence. Not only was the 
permeability sequence different, we also found that the 
Ici.ca of mTMEM16F-expressing cells had a significantly 
higher Pbi/Pci than the Ici.ca in mock-expressing cells. 

Anion selectivity of mTMEM16F mutants 

A previous study on TMEM16A mutants (Yang etal., 
2008) indicated that the pore region lies within the highly 
conserved extracellular domain between TM5 and TM6. 
Mutation of positively charged amino acids R621, K645, 
and K668 to the negatively charged glutamic acid within 
this region resulted in changes in the kinetics of anion 
currents and a reduction in the anion/cation permeabil- 
ity ratio. To investigate the function of this region in 
TMEM16F, we introduced mutations in R592, K616, and 
R636 (Fig. 6 E) by replacing them with glutamic acid (E) 
and constructed a control mutation I342A which we an- 
ticipated would not affect the conductive properties of 
the putative channel pore. We introduced the R592, K616, 
and R636 mutations into both GFP coexpression vec- 
tors and the C-terminal GFP fusions. Fig. 6 A shows the 
peak Ici.ca activation of HEK293 cells transfected with 
mock, mTMEM16F, or the single amino acid mutants of 
mTMEM16F. Compared with the maximal whole-cell Ca 2+ - 
activated (250 uM) CP current in mTMEMl 6F-WT, the 
mutation R592E reduced the current strongly and was no 
longer significantly different from currents in mock trans- 
fectants. In contrast, the maximal current was not signifi- 
cantly affected in the I342A, K616E, and R636E mutants 
(Fig. 6 A) . Fig. 6 (B-D) and Table 1 show the relative 
anion permeabilities recorded in whole-cell mode after 
stimulation with 250 uM [Ca 2+ ] ; for I342A-, K616E-, and 
R636E-expressing HEK293 cells, respectively. The Psr/Pci 
of the control mutation I342A was not significantly differ- 
ent from mTMEMl 6F-WT (Table 1). However, the obser- 
vation that in mTMEMl 6F-WT P Br was significantly higher 
than P C i was not reproduced with the I342A control mu- 
tation, perhaps as a result of the variability in the data. 
Although the permeability sequence was not changed in 
K616E, the relative permeability of iodide was lower 
when compared with mTMEMl 6F-WT (Table 1). In 
R636E, a marked reduction was found in the relative per- 
meability to thiocyanate, iodide, and bromide (Fig. 6, 
C and D; and Table 1) when compared with those found 
in mTMEM 1 6F-WT (Table 1). 

Subcellular localization of mTMEMl 6F mutants 

The subcellular localization of the R636E, K616E, and 
R592E single amino acid mutants are shown in Fig. 6 F. 
Prominent staining is seen in the plasma membrane as well 
as in intracellular compartments of the mTMEM16F-WT 
as well as mutant-expressing cells (Figs. 2 and 6 F) , suggest- 
ing that the loss of current in the mutants is not caused by 
trafficking problems. To investigate the findings from the 
immunofluorescence imaging further, we performed a 
biotinylation assay in which plasma membrane proteins 
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Figure 5. Anion selectivity of the 
mTMEM16F-induced current. (A and 
B) Representative anion selectivity cur- 
rent traces recorded in whole-cell con- 
figuration after stimulation with 250 uM 
[Ca 2+ ]j of mock (A)- and mTMEM16F- 
transfected (B) HEK293 cells. (C and 
D) Mean relative anion permeabilities 
(Px/Pci) for mock (Q- and mTMEM16F- 
transfected (D) HEK293 cells. P x /P a are 
presented as a mean of seven to nine cells 
and are tested by one way ANOVA; "*" 
indicates significant (P < 0.05) differ- 
ence from PCI, and "f" indicates signi- 
ficant difference (P < 0.05) between the 
relevant anion permeabilities when 
tested with paired t test. Only experi- 
ments in which the anion permeabilities 
were paired were included in the statisti- 
cal test. 



were first isolated and TMEM16F then visualized by West- 
ern blotting. The membrane proteins were isolated from 
HEK293 cells transfected with plasmids containing the 
respective TMEM16F-GFP constructs (WT and the four 
mutants I342A, R592E, K616E, and R636E). As seen in 
Fig. 6 G, mTMEM16F-WT along with all of the mutants 
localizes to the plasma membrane, consistent with the re- 
sults from the immunofluorescence imaging. 

The Na + current component of membrane currents 
in TMEM16F-expressing cells 

The protocols applied above allowed for a quantitative 
analysis of anion selectivity ratios that were independent of 
the plasma membrane's Na + permeability, however large it 
might have been. Here, we analyze absolute Na + and CP 
permeabilities of HEK293 cells expressing WT TMEM16F 
and the various mutants of the channel. Fig. 7 A shows 
instantaneous and (near) steady-state current-voltage rela- 
tionships of a TMEM16F-WT-expressing HEK293. With 
NMDG + in bath, the instantaneous current, I Inst , could be 
fitted by the GHK CP current equation 



J Ins, -p 



F 2 V m [Cr\Exp(-FV m /RT) -[GI- 



RT 



Exp{-FV m / RT)-l 



(3) 



In Eq. 3, [CI ]; = 44.7 mM and [CT] 0 = 155 mM; i.e., the 
CP equilibrium potential is E a = —31 mV. With the CP 
permeability as free parameter, we obtained Pci = 12.9 
(±0.6) x 1(P 14 cm 3 /s/pF Replacement of NMDG + with 
Na + resulted in rightward displacement of the I Inst /V m re- 
lationship as expected if the membrane current also con- 
tains an Na + component. The full line of Fig. 7 A (left) 
indicates the nonlinear fit of Eq. 4, which expresses the 
instantaneous currents as the sum of two GHK-current 
components carried by CP and Na + , respectively: 



J Ins, _ p 



RT 



Exp(FV m /RT)-l 



(4) 



F 2 V m [Cr\Exp{-FV m /RT)-\Cr\ 
RT Exp(-FV m / RT)-\ 

F 2 V m [Na + \Exp{FV m / RT) - [Na + ] 0 



TABLE 1 

Anion permeability ratios relative to P a of mock-, mTMEM16F-, I342A-, K616E-, and R636E-transfected HEK293 cells 



Px/Pci 


Mock 


TMEM16F 


I342A 


K616E 


R636E 


Pscn/Pci 


4.20+1.03 (NS) 


3.46 ± 0.80 


2.62 ± 0.38 (NS) 


2.19 ±0.23 (NS) 


1.64 ±0.37* 


Pi/Pq 


1.47 ±0.16 (NS) 


1.82 ±0.20 


1.85 ±0.22 (NS) 


1.35 ±0.11" 


0.86 ± 0.09* 


P Br /Pcl 


0.96 ±0.11* 


1.32 ±0.06 


1.32 ±0.21 (NS) 


1.09 ±0.08 (NS) 


0.86 ±0.11* 


Pa, p /Pci 


0.57 ± 0.08 (NS) 


0.51 ± 0.06 


0.53 ±0.11 (NS) 


0.70 ± 0.06 (NS) 


0.50 ± 0.06 (NS) 



The ratios were calculated by Eq. 2. 

Significant differences between relative anion permeabilities of mTMEM16F-WT and mock or mutant mTMEM16F-expressing cells are indicated 
(P < 0.05), when tested with one-way ANOVA and using a Student-Newman-Keuls post hoc test. 
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Figure 6. Effect of single amino acid mutations on whole-cell current and anion selectivity. (A) Mean peak Ic a ,ci current density when 
activated by 250 uM free [Ca 2+ ]i of mock, mTMEM16F, and single amino acid mutants in the putative pore site of mTMEM16F (the 
control mutation I342A and the pore region mutants R592E, K616E, and R636E) measured at 70 mV and — 1 10 mV, respectively. Data 
are presented as mean ± SE; you can read the amount of replicates inside the 1 10-mV bar graphs; "*" and "f" indicate significant (P < 
0.05, when tested with one-way ANOVA) difference from mock and mTMEM 1 6F, respectively. (B-D) Anion permeabilities recorded in 
whole-cell mode after simulation with 250 pM free [Ca 2+ ]; for I342A (B)-, K616E (C)-, and R636E-transfected (D) HEK293 cells. Data 
are presented as mean ± SE; the number of replicates are indicated in the figures; "*" indicates significant (P < 0.05) difference from 
PCI, and "f" indicates significant difference (P < 0.05) between the relevant anion permeabilities, when tested with paired itest. Only 
experiments in which the anion permeabilities were paired were included in the statistical test. (E) Putative topology of the TMEM16F 
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With [Na + ] ; = 15.3 mM, [Na + ] Q = 150 mM, and the two ion 
permeabilities being free parameters, we obtained Pci = 
31.7 (±1.0) x 10~ 14 cm 3 /s/pF and P Na = 9.7 (±0.5) x 10~ 14 
cm 3 /s/pF. The Fig. 7 A (right) shows the corresponding 
outwardly rectifying (near) steady-state current-voltage 
relationships. I StSt with NMDG-C1 in bath intersected the 
x axis at V rev ~ E a . However, with NaCl in bath and E Na ~ 
57 mV, the reversal potential was displaced to the right 
of E C i, in casu V rev ~ — 12 mV, confirming that the con- 
ductance of the plasma membrane is governed also by a 
significant Na + permeability (Fig. 7 A, inset) . 

Similar acceptable nonlinear fits of Eqs. 3 and 4 were 
obtained with all other mutants. This is exemplified in 
Fig. 7 B with HEK293 cells expressing the TMEM16F- 
I342A, mutated at a site outside the pore region, and in 
Fig. 7 C with TMEM16F-R636E, mutated in the pore re- 
gion which significandy altered the anion selectivity ra- 
tios. Table 2 collects the results of the analysis of all I Inst /V m 
relationships generated by the above protocols. There 
was no systematic variation of PNa/Pci between the groups, 
indicating that the ratio is not affected by the mutations 
introduced here. As shown Fig. 8, cells with large Pci also 
exhibited large Pn 3 . The cross correlation coefficient ob- 
tained by a least squares fitting to all data was R~ = 0.903, 
with a linear relationship between the two permeabilities 
given by P Na = 0.334 x P cl (±SD = 0.028; n = 15). 

Comparing whole-cell CP permeabilities obtained with 
NMDG + and Na + in bath showed that in 13 of the 15 cells 
studied, Pci increased significantly by replacing NMDG + 
with Na + . This result was independent of the expressed 
channel construct. The mean ratio (±SEM) of Pci with 
NMDG + outside and P C i with Na + outside was 0.64 (±0.07; 
Table 2). 

DISCUSSION 

The current carried by TMEM16F 

Mutations in genes of the subfamily among the human 
TMEM1 6 proteins (TMEM1 6E and TMEM 1 6F) have been 
associated with gnathodiaphyseal dystrophy (Katoh and 
Katoh, 2004; Tsutsumi etal., 2005) and Scott syndrome 
(Suzuki etal., 2010; Castoldi etal., 2011), respectively, 
and TMEM16F has been shown to be necessary for the 
Ca 2+ -dependent phospholipid scramblase activity (Suzuki 
et al., 2010; Yang et al., 2012). The literature about the 
putative function of TMEM16F as an ion channel is con- 
tradictive. Caputo et al. (2008) found that siRNA knock- 
down of TMEM16F did not affect the anion current in 
their halide-sensitive YFP assay, and Duran et al. (2012) 



recently published that TMEM16F is an intracellular 
protein. In contrast, although not characterized in detail, 
other studies have suggested that TMEM16F produces 
Ca 2+ -sensitive currents in FRT and HeLa cells (Schreiber 
etal., 2010; Kunzelmann etal., 2012). To complete the 
confusion, it has been published that TMEM16F forms a 
Ca 2+ -activated nonselective cation channel (Yang etal., 
2012) . In the present study, we characterized the electro- 
physiological properties of the TMEM16F-associated an- 
ion currents in WT TMEM16F and in constructs mutated 
in the putative pore region. Moreover, we investigated 
the localization of TMEM16F and the relation between 
P Na and Pci after activation of TMEM16F. 

TMEM16F localizes to the surface membrane 
Although we noted relatively strong mTMEM16F stain- 
ing of intracellular organelles, we also observed dis- 
tinct mTMEM16F localization to the plasma membrane 
in HEK293 cells (Figs. 2 and 4). This is consistent with 
previous localization studies of TMEM 16 family mem- 
bers including TMEM16A, TMEM16B, TMEM16E, and 
TMEM16G (Bera etal., 2004; West etal., 2004; Rock 
et al., 2008; Schroeder et al., 2008; Yang et al., 2008). 

Ca 2+ simultaneously activates with delay voltage- 
and time-dependent anion and cation currents 
in mTMEM16F-expressing cells 

We have shown that over expression of mTMEM16F in 
HEK293 cells results in [Ca 2+ ]i-regulated CP currents 
and a significant cation conductance. Furthermore, the 
observed anion activity shows significandy different char- 
acteristics to what has been reported for both TMEM16A 
and TMEM16B as well as native CaCCs in various tissues 
(Caputo etal., 2008; Schroeder etal., 2008; Yang etal., 
2008; Pifferi et al., 2009; Stohr et al., 2009; Schreiber 
et al., 2010). The current was not activated immediately 
after establishment of whole-cell configuration with high 
[Ca 2+ ] in the pipette solution. However, after a delay of sev- 
eral minutes, a transiendy activated and strongly outwardly 
rectifying current that peaked after a mean of ^9 min 
was observed. In the study by Schreiber et al. (2010) , a 
slowly activating current was reported for TMEM16F- 
expressing FRT cells, which has littie resemblance with 
our observations, perhaps because their observation pe- 
riod was no more than 8 min. However, in a recent study 
by Shimizu etal. (2013), treatment with ionomycin in- 
duced a delayed and slowly activating current, which is 
similar to our findings. In our study, mTMEM16F was 
characterized by a strong outwardly rectifying current in 
asymmetrical solutions, a time-dependent activation at 



channel with the mutated amino acids' approximate positions. (F) Immunofluorescence images of localization of TMEM16F mutants 
expressed in HEK293 cells. The images were taken using an epifluorescence microscope. GFP was enhanced by anti-GFP and Alexa 
Fluor 488 (green). (G) Presence of mTMEM16F-GFP in the plasma membrane fraction and total cell lysate of HEK293 cells expressing 
either GFP-tagged WT mTMEM16F or one of the mutants I342A, R592E, K616E, and R636E. Plasma membrane proteins were isolated 
using a biotinylation assay. [5-Actin and NPTII were used as controls for loading and transfection efficiency, respectively. 
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Figure 7. Instantaneous (I Imt ) and near steady-state (I StSt ) current- 
voltage relationships of HEK293 cells in whole-cell patch-clamp 
configuration. Green symbols show I-V relationships obtained 
with NMDG-CI in bath. Blue symbols show I-V relationships ob- 
tained with NaCI in bath. I Inst was recorded immediately after the 
voltage step from the holding potential, Vhoid = 0 mV, to the new 
V m . I St " St was recorded toward the end of the 1-s voltage pulse. 
Thus, rather than estimating the PNa/Pci of each cell from a rever- 
sal potential obtained by interpolation, we could take advantage 
of I Inst -V m data points covering a broader range about V rev . Cur- 
rent-voltage relationships are shown for the following (all perme- 
abilities in units of 10~ 14 cm 3 /s/pF) . (A) HEK293 cells expressing 
TMEM16F-WT (with NMDG-CI in bath, P cl = 12.9 ± 0.6; with NaCI 
in bath, P cl = 31.7 ± 1.0 and P Na = 9.7 ± 0.5). (B) HEK293 cells 
expressing the neutral mutant TMEM16F-I342A (with NMDG-CI 
in bath, P cl = 19.8 ± 0.8; with NaCI in bath, P cl = 93.8 ± 4.2 and 
?Na= 13.3 ±3.6). (C) HEK293 cells expressing the anion sequence 
modifying TMEM 1 6F-R63 6E mutant (with NMDG-CI in bath, 
P C1 = 16.9 ± 1.5; with NaCI in bath, P cl = 22.9 ± 1.7 and P Na = 8.0 ± 
1.4). It is noted that for all three I M /V m relationships obtained 



depolarizing potentials (positive to 40 mV), and a time- 
dependent deactivation at hyperpolarizing potentials 
(Vhoid = ~~ 20 mV), which resulted in large steady-state 
outward currents and relatively very small steady-state 
inward currents, similar to what was found for TMEM 1 6B 
(Pifferi etal., 2009; Stohr etal, 2009). The very slow 
activation of Ca 2+ is not yet explained. We considered 
that it was caused by delayed trafficking, for example if 
recruitment of TMEM16F to the plasma membrane de- 
pends on a slow Ca 2+ activation of a second messenger. 
This hypothesis was tested; however, under the condi- 
tions of our experiments, we did not find support for 
this hypothesis (Fig. 4 E). Shimizu etal. (2013) investi- 
gated more specifically whether CAMKII has a role in 
this delay, but they failed to see any effect using a CamKII 
inhibitor (KN93) , and the very slow activation is yet to 
be explained. 

Ca 2+ activation kinetics of mTMEM16F 

We found that the delayed Ia,ca activations in mTMEM16F- 
expressing cells depend on [Ca 2+ ]; in the micromolar 
range (Fig. 4, A-D). The EC 50 for Ca 2+ is —76 pM, and 
when measured at depolarizing potentials, the Hill coef- 
ficient was calculated to 2, in line with Yang et al. (2012) , 
suggesting cooperative Ca 2+ activation (Fig. 4 C). In 
comparison, an EC50 of 0.4 pM and 3.3 pM with an ra H 
of 2 and >2 has been reported for TMEM16A and 
TMEM16B, respectively (Yang etal., 2008; Pifferi etal., 
2009; Scudieri et al., 2012) . For TMEM16F, it was recendy 
reported that the EC50 value was 3.4 pM (Yang etal., 
2012) and 9.6 pM (Shimizu etal., 2013). Both groups 
have used EGTA-buffered solutions and calculated the 
free Ca 2+ concentration using Cabuf software. Using the 
same software, we find an EC50 value of 12 pM Ca 2+ . Be- 
cause the Ca 2+ concentrations of use are well outside the 
dynamic range for the EGTA buffers (Dweck et al., 2005) , 
we have measured the actual Ca 2+ concentration of the 
pipette solution (see Materials and methods), finding a 
significandy higher free Ca 2+ concentration, shifting the 
dose-response curve rightwards. Thus, die difference in 
Ca 2+ dependence between the three groups most likely 
reflects methodical difference in [Ca 2+ ]f ree estimation 
rather than significandy different affinities. However, the 
Ca 2+ dependence is dramatically different from the results 
for TMEM16A (compare also with Fig. 3 C) . 

The slope of the I/V relationship of the strongly out- 
wardly rectifying TMEM16F CP current did not approach 
linearity at increasing intracellular Ca 2+ concentrations 
(Fig. 4 A) , as seen for endogenous CaCC (Kuruma and 
Hartzell, 2000), TMEM16A (Schroeder etal., 2008), 
and TMEM16B (Pifferi etal., 2009) currents. Contrary 



with NMDG-CI in bath that V rev » E a , whereas the reversal poten- 
tials with NaCI in bath are significantly displaced to the right of 
E C i, which is emphasized by the insets of expanded axis scales. 
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TABLE 2 

Cr and Na* permeabilities obtained from GHKfits to instantaneous current-voltage relationships of 15 HEK293 cells 
Cell type NMDG-C1 in bath NaCl in bath P N a/Pci PcTNMDG-a/Pcf Naa 





Pa 


±SD 


Pa 


±SD 


l\,i 


±SD 








ltT 14 cm 3 /s/pF 


10-" cm 3 /s/pF 


Itr 14 cni'/s/pF 


1CT 14 cm 3 /s/pF 


ltT 14 cm 3 /s/pF 


KF 14 cm 3 /s/pF 






WT 


30.4 


1.9 


29.6 


1.2 


7.1 


0.6 


0.24 


1.03 


WT 


12.9 


0.6 


31.7 


1.0 


9.7 


0.5 


0.31 


0.41 


I342A 


26.8 


0.8 


55.6 


2.0 


14.7 


1.7 


0.26 


0.48 


I342A 


19.8 


0.8 


93.8 


4.2 


13.3 


3.6 


0.14 


0.21 


I342A 


138.5 


15.0 


230.2 


1.6 


77.9 


1.4 


0.34 


0.60 


I342A 


87.2 


6.9 


127.6 


2.4 


67.7 


1 .9 


0.53 


0.68 


K616E 


204.1 


16.0 


176.8 


12.7 


45.8 


10.6 


0.26 


1.15 


K616E 


45.7 


4.7 


137.7 


7.7 


30.6 


6.4 


0.22 


0.33 


K616E 


12.8 


1.0 


23.4 


2.5 


12.7 


2.1 


0.54 


0.55 


R592E 


10.6 


3.4 


12.7 


2.3 


9.7 


2.0 


0.76 


0.83 


R592E 


42.4 


1.3 


72.7 


3.5 


11.2 


3.0 


0.15 


0.58 


R636E 


16.2 


1.7 


27.3 


1.8 


3.1 


1.5 


0.11 


0.59 


R636E 


90.0 


2.4 


130.4 


1.2 


69.5 


1.0 


0.53 


0.69 


R636E 


16.9 


1.5 


22.9 


1.7 


S.O 


1.4 


0.35 


0.74 


R636E 


150.0 


8.5 


232.0 


4.6 


79.9 


3.8 


0.34 


0.65 



As explained in Eq. 4 (compare with Fig. 7, A-C). The error of each estimate is indicated by the SD of the fit given by the Origin fitting routine (OriginLab). 

TMEM16F-expressing HEK293 cells conduct a signifi- 
cant Na + current and that cells expressing the mutated 
TMEM16F exhibited a significant Na + permeability as 
well. Moreover, within the range of Pci from ^1 x 10 
to 25 x 10 cm/s, the two permeabilities were corre- 
lated according to P Na = 0.334 (±0.028) x P cl (Fig. 8). 
In the face of the anion selectivity series obtained, it is 
unlikely that Na + shares the TM pore with CP. Alterna- 
tive interpretations would be that trafficking of Na + and 
CP channels to the surface membrane is coupled or, 
perhaps, that the TMEM16F molecule expresses paral- 
lel CP and Na + selective pores. In connection to this, it 
is interesting that TMEM16F has been described both 
as a Ca 2+ -activated nonselective cation channel with 
PNa/Pci = 6.8 (Yang etal., 2012) and as a Ca 2+ -activated 
selective anion channel (Shimizu et al., 2013) . Although 
Yang et al. (2012) was unable to observe any shift in rever- 
sal potential upon shifting from CI to the larger anion 
Mes, Shimizu etal. (2013) did not find any change 
when shifting Na + cations to the bigger NMDG ions. In 
contrast, we observed a significant shift in E rev , both 
when substituting Na + with NMDG and when substitut- 
ing CP with aspartate. The background of this differ- 
ence is not clear, but an important point is that mutations 
in TMEM16F have been observed to change the perme- 
ability profile for anions (K616E and R636E; this study) 
and the measured cation current (Q559K; Yang et al., 
2012), indicating that the observed currents are carried 
by TMEM16F. However, it is unclear whether the decrease 
in PNa/Pci observed for the Q559K mutation (Yang 
et al., 2012) actually is a decrease in P Na or potentially 
could relate to an increased CI permeability. Thus, it 
is so far unclear whether the selectivity filter for cat- 
ions is to be found in the TM5-6 region as seen for the 



to what has been reported for TMEM16A (Schroeder 
etal., 2008; Yang etal., 2008) and TMEM16B (Pifferi 
et al., 2009) , the TMEM16F-generated current shows time- 
dependent activation at depolarizing potentials even at a 
saturating free [Ca 2+ ] j, (Fig. 4, A and D) . This would indi- 
cate that full activation of the mTMEM16F-generated 
current is not governed solely by [Ca 2+ ];. 

The Na + ion permeability of TMEM16F-expressing cells 
Surprisingly, replacing extracellular NMDG + with Na + 
stimulated the CP permeability in a large majority of 
studied TMEM16F constructs (Table 2). The mechanism 
of the interaction between extracellular Na + and the 
CP current is not known. Our study also showed that 




P Ci (10" 7 cm/s) 

Figure 8. Relationship between P a and P Na of HEK293 cells ex- 
pressing WT and engineered TMEM16F constructs. The perme- 
abilities were obtained as described in the legend of Fig. 7. For 
expressing permeability in the conventional units of cm/s, the 
values given in Table 2 have been recalculated using the empiri- 
cal approximation that 10 6 pF corresponds to a membrane area 
of 1 cm 2 . The linear fit to all data is given by P Na = 0.334 x P cl 
(±SD = 0.028) , R 2 = 0.903, and n = 15. 
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anion selectivity filter (Fig. 6; Yang et al., 2008) . This will 
have to be further tested in the future, preferably by sin- 
gle channel analysis. 

Anion selectivity of mTMEM16F-WT and mTMEM16F 
mutants: TMEM16F is a functional anion channel 

The anion substitution experiments showed that the 
Ici.ca of mTMEM16F follows the Eisenman type 1 selectiv- 
ity sequence, P S cn > Pi > Pbf > Pci > P^p; indicating weak 
interactions of the permeating ion with the channel 
pore, which has been reported also for CaCCs and volume- 
regulated anion channels (Nilius and Droogmans, 2003), 
TMEM16A (Yang et al., 2008), and TMEM16B (Pifferi 
etal., 2009). 

To investigate whether the TMEM16F protein constitutes 
a channel for halide ion currents, we performed mutations 
of amino acids in the putative pore region between TM5 
and TM6. Within this region, mutations in TMEM16A 
of the positively charged amino acids R621, K645, and 
K668 to the negatively charged glutamic acid resulted in 
diminished selectivity ratios relative to CP for several an- 
ions (Yang et al., 2008) . We introduced the following single 
amino acid mutations, R592E, K616E, and R636E, replac- 
ing the positively charged arginine-592, lysine-616, and 
arginine-636 with the negatively charged glutamic acid. Al- 
though these mutants all trafficked to the membrane, the 
mutant R592E reduced the current strongly, whereas the 
maximal current was not affected in the K616E and R636 
mutants (Fig. 6 A). Notably, we found a reduced Pi/Pa in 
K616E when compared with mTMEM16F-WT (Table 1) 
and an anion selectivity sequence changed to Pscn = Pi = 
P Br = Pa > P% P in the R636E mutant (Fig. 6 C). The above 
results suggest that these residues are involved in the anion 
selectivity of the pore. In contrast, the current density and 
anion selectivity sequence was unaffected in the control 
mutant I342A, in which the amino acid substitution was 
located outside the putative pore region (Fig. 6, B and E) . 
Collectively, these findings provide the evidence that 
TMEM16F is a pore-forming subunit of an anion channel. 

Conclusion 

Mouse TMEM16F heterologously expressed in HEK293 
cells localizes to the plasma membrane and induces 
an Ici.ca with properties different from canonical Ici.ca 
and from Ca 2+ -activated anion currents associated with 
TMEM16A and TMEM16B. Furthermore, a significant 
Na + current was activated, and the two permeabilities were 
correlated according to Pn 3 = 0.3 Pci- 
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